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O ' Abstract 

m 

0> ■ We study the phase diagram of statistical systems of closed and open in- 

terfaces built on a cubic lattice. Interacting closed interfaces can be written 
as Ising models, while open surfaces as Z{2) gauge systems. When the open 
surfaces reduce to closed interfaces with few defects, also the gauge model 
I can be written as an Ising spin model. We apply the lower bound renormal- 

ization group (LBRG) transformation introduced by Kadanoff (Phys. Rev. 
O ■ Lett. 34; 1005 (1975)) to study the Ising models describing closed and open 

surfaces with few defects. In particular, we have studied the Ising-like tran- 
sition of self-avoiding surfaces between the random-isotropic phase and the 



X 



5_j ■ phase with broken global symmetry at varying values of the mean curva- 



ture. Our results are compared with previous numerical work. The limits of 
the LBRG transformation in describing regions of the phase diagram where 
not ferromagnetic ground-states are relevant are also discussed. 



PACS number: 68.10.-m (Fluid surfaces and fluid-fluid interfaces); 05.50.+q (Lat- 
tice theory and statistics; Ising problems); 64.60.Ak (Renormalizat ion-group, frac- 
tal, and percolation studies of phase transitions); 11. 15. Ha (Lattice gauge theo- 
ries) ; 



§1. Introduction. 

In this paper we will apply a renormalization group transformation to study 
the phase diagram of interface models built on a cubic lattice. Fluid interfaces 
in 3D statistical systems are subject of much current research They provide 
useful descriptions of experimental systems such as mixtures of oil, water and 
surfactant, or aqueous solutions of surfactant 0. In ternary mixtures the sur- 
factant forms monolayered interfaces between oil and water; in aqueous solutions 
bilayered membranes are typical constituents of biological cells. The properties 
of these systems at low surfactant concentrations are relevant for both practical 
and theoretical reasons. For example, in ternary mixtures, a middle phase 
coexisting with oil-rich and water-rich phases is considered very appealing for ap- 
plications 1^, due to the very low surface tension values between the coexisting 
phases. From a theoretical point of view, dilute interfaces can be seen as ex- 
perimental realizations of random surface models where self-avoidness is the only 
relevant interaction 

The typical lack of topological constraints on the physical configurations sug- 
gests the use of lattice models to describe ensembles of fluid surfaces. First con- 
sider the case of closed interfaces without defects such as holes or seams. Closed 
interfaces can be described in Ising models as the boundaries separating domains 
of opposite spins, which, in the identification with ternary mixtures, can represent 
oil and water. The interfaces are built on the dual lattice, and, for a given spin 
configuration {<7i}, have a total area S = J2<ij>{^ — c"jcrj)/2, where the sum is 
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over all nearest-neighbour pairs in the original lattice. Other surface energies can 
be considered by introducing further spin interactions. Surfaces where curvature 
and intersections are also weighted can be represented by a generalization of 
the Ising model defined by the hamiltonian [7-9] 

H = JiY^ aidj + J2 ^i^j + ^4 X! (yi(yj(ykcri (1) 

<ij> W 

I k 

where the three sums are respectively over nearest, next-to-the- nearest neigh- 
bours and plaquettes of a cubic lattice. Here and in the following of the paper, 
our definitions of hamiltonian will always include tha factor — /3. The parameters 
J\-,Ji-,J^ can be expressed in terms of the surface parameters (isificii^L H] rep- 
resenting respectively the energy cost for an elementary area (one plaquette on 
the dual lattice), for two plaquettes at right angle, and for four plaquettes with a 
common bond [|T^. A positive favours fiat configurations; it corresponds to a 
mean curvature energy, which has been proved to be an useful phenomenological 



parameter for describing fluid interfaces |TT[]. The term proportional to /3l can 
mimic the self-avoidness interaction in the limit /S^ — 00, when surfaces touching 
each other along some contour are forbidden. The relations between spin and 
surface parameters are: 

Ji- ^ + /^c, ^2 = -Y-^, ^4 = -^ + ^. (2) 

The phase diagram of the model (|I|) has been studied by mean-field and numerical 
simulations in [§, |^ . It exhibits many properties relevant for real systems, as also 



discussed in [|I2 



If one wishes to consider the effects of defects in fluid interfaces, ensembles of 
open surfaces have to be introduced. A simple lattice realization of open surfaces 



T3| is given by the self-dual Z{2) gauge model [T^. Here 2- values variables {Uij} 
are defined on the bonds of a cubic lattice. One says that the plaquette dual to the 
bond < ij > is occupied by some surface if Uij = —1; it is not occupied if Uij = 1. 
Therefore a given {Uij} configuration corresponds on the dual lattice to a surface 
configuration with area S = J2<ij>{^ ~ Uij)/2. A bond on the dual lattice can 
be said to belong to some defect if an odd number of the dual plaquettes sharing 
that bond is occupied by some surface. Defects defined in this way can be counted 
by considering the product of Uij over the bonds of each plaquette in the original 
lattice. It is easy to recognize that the total length of defects will correspond 
to the quantity V = J^i^i i.^ ~ UijUjkUkiUii)/2 |13[. Therefore the self-dual Z{2) 
gauge model with hamiltonian given by 



(3) 



H = Pm E + /5g E Ui,U,kUkiUu 
<ij> 

describes open surfaces where area and defects are both weighted. Self-duahty 



here means that the model is symmetric with respect to the transformations 



^ (^G = -^Intanh/^M, 



/3g ^ /^M = -- In tanh pc- (4) 



In the parametrization (|^) a large value of Pm favours configurations with small 
area, while a large value of jSc inhibits defects. 

The phase diagram of the self-dual Z{2) gauge model has been first analyzed 
1^ , pTSU ; it has been studied by Monte Carlo simulations in At small 
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values of /3m - and analogously, by self-duality, at large values of (3g - it can be 



shown |T5|, [13] that the model @ can be expanded as an Ising spin model with 
an increasing number of interactions. For example, at the second order of the 
expansion at small (3m the model (0) can be written as 

H = aiaj + J2 ^ CTjCTj + J4 51 o-iajakcri + Je X! (^i(^j(^k(^i (5) 

<ij> irf' cor 

I k 

where the interactions are between nearest neighbours, next-to-the-nearest neigh- 
bours, the 4 spins of a plaquette and the 4 spins of a corner (see Table 1). The 
coupling constants J\,j2-,Jii, and jTe can be expressed in terms of the constants 
/3g and /?m as follows 

= (3- 4(tanh /?m)^[3 cosh 2/3 sinh 2(3 + (cosh 2/3)^ sinh 2(3] 
J2 = 2(tanh /5m)'' [(sinh 2/3)2 + cosh 2/3(sinh 2/3)2] 

(6) 

J4 = 2(tanh/3M)^(sinh 2/3)2 
Je = -i(tanh/?M)''(sinh2/3)3 



where 



(3= -- Intanh [/3g + (tanh /3m)^] • (7) 
2 



The spin representation has the advantage that it can be more easily studied . 



In this paper we will apply the so called lower-bound renormalization group 



(LBRG) transformation first proposed by Kadanoff |T9[ to study the phase dia- 
gram of the models (0) and (^ as given in the approximation (|^). The LBRG 
transformation can be conveniently applied to cases where all the interaction is 
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in an elementary cell of the lattice, as it is in the models ([^) and @. The conve- 
nience is appreciable especially in D = 3 where other RG transformations would 
be much more dispendious from a computational point of view. 

The LBRG transformation produces a lower bound to the free-energy which 
can be maximized by conveniently fixing a variational parameter. Its application 
to various models generally gives very accurate estimates of critical exponents 



For example, in the 2D Ising model it predicts the inverse critical temperature 
Peril = 0.458 iPfrT'^"'' = 0.4407), and the exponent u = 0.999 (i/On^ager ^ 



The drawback of the LBRG transformation is that it preserves the nature of the 
ground states only in the ferromagnetic region, so that it can be reliably applied 
in a limited region of the phase diagram. 

We will describe the LBRG transformation in §2. In §3 we will show the results 
obtained by applying the LBRG transformation to the models (|1|) and (0,^. In 
particular, in model (|I|), the self-avoidness limit is examined for different values 
of the curvature Pc- A discussion of our results will follow in §4. 



§2. The LBRG Transformation. 

In this section we will briefly describe the LBRG transformation. Further 
details can be found in Real space RG transformations can be generally 
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written as 



e 



Here J denotes a set of coupling constants, o = {cri, ...(Jn} a spin configuration 
and H{a, J') is the hamiltonian to be studied; the weight function V{(t'; cr) defines 
the renormahzed hamiltonian H'{a', J'') with new spin variables a'^, {N' < 

N) and coupling constants J''. The relation J^a' ^(c^'i o") = 1 ensures that the 



total free-energy is unchanged. In the LBRG transformation the spin o"^ are 
defined on the cells like those marked by a cross in Fig.l; P(cr', a) is chosen as the 
product over the marked cells of the functions 

2cosh[p(cri,i + ... + ai^8)\ 

with p a real parameter and CTj^i, ...cTj^s the original spins at the vertices of the 
cube i. 

If the original hamiltonian can be written as H{a,J') = J2 H{(J, JC), where 
the sum is over the elementary cubes of the lattice and /C is the set of couplings 



normalized to a single cell a convenient moving of interactions and factors of 
d^) will give a new hamiltonian with all the interaction still in a single cell. The 
renormahzed cell hamiltonian H'{(t', /C') is given by 

r/v// / / WM exp[p{a[(ri + ... + a'gas) + 8H{ai,...,as;IC)] 
exp|ff(.,.....,;C)]=^£^ . 

(10) 

Since we are interested in studying the phase diagrams of the Ising models (|l|) 
and (I) where only even interactions appear, it will be sufficient to consider the 



transformation of the 14 even couplings (see Table I) which can be defined on a 
cell of a 3D cubic lattice. One of these couplings is a pure constant, we denote it 
by /Co; the others are denoted by /Cj z = 1, 13. 

After some algebra one gets from (llOl) the recursion laws 



\ IC'q= 8/Co + $o(p;/Ci,...,/Ci3) . . 

\ 1C[= <l>,(p;/Ci,.../Ci3) V^ = l,...,13 

where $0) •••) '^'is are analytic functions. The critical properties of the system can 

be then related to the behaviour of the recursion laws close to their fixed points. 

The variational nature of the interaction-moving operation was first observed 



by Kadanoff ||2^. A lower bound f*{p) to the free-energy per site can be calculated 

by 

f*(p) = - hm ^ ; (12) 



where /Cg is the value of /Cq after n applications of the LBRG transformation [O 



Following the prescription of [|I9[, the parameter p will be fixed by maximizing 
the function f*{p) starting the iterations from the fixed point hamiltonian with 
/Co = 0. 



§3. Results. 

Closed interfaces - model The LBRG transformation is here applied to 
calculate the ferromagnetic-paramagnetic (F-P) transition surface in the space 
JxiJiiJi^- For completeness, results concerning other transitions, related to not 
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ferromagnetic ordering, will be also given. These results have to be considered 
with cautions since the LBRG transformation, as defined in §2, does not take 
correctly in account the structure of not ferromagnetic ground states. 

The value of p maximizing the critical fixed point free-energy on the F-P 
surface is p* = 0.40354. In Table 1 fixed points related to the F,P and AF 
(antiferromagnetic) phases are reported for the value p = p*- 

The fixed points (f), (p) and (c) are respectively the low-temperature ferro- 
magnetic, the high temperature and the F-P critical fixed points. The LBRG 
transformation has been already applied for calculating the exponents of the 3D 
Ising model in JTol, where the optimal value found for p* is p = 0.40343. We do 



not understand the reasons of the discrepancy with our result. The critical fixed 
point at p = 0.40343 is reported in the caption of Table 1. At p = p* the values 
of the inverse Ising critical temperature and of the exponent u are respectively 
Pcrit = 0.23925 and u = 0.6288. The corresponding values at p = 0.40343 are 
f3j.^n = 0.23923 and u = 0.6290; the best estimates are (3^^ = 0.22165 and 
u = 0.6289 ± 0.0008. The fixed point (c) is symmetric in the sense that all the 



2-spin, the 4-spin, etc. interactions are equal. This symmetry was assumed in [p!9 
while here we consider recursions in the whole space of couplings. This situation 
can be compared with the results obtained by applying the LBRG transformation 



to the 2D Ising model [21]. In D = 2 the symmetric fixed point has two relevant 
eigenvalues with an eigenvector pointing outside the symmetric subspace on the 
critical surface. Therefore in D = 2, differently from the 3D case, the symmetric 
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critical fixed point, which is found to maximize the free-energy, cannot be reached 



starting from non symmetric interactions |21[ . 

The fixed points (af) and (ac) are the antiferromagnetic counterparts of the 
fixed points (f) and (c). The fixed point (ac) is on the transition surface between 
the AF and the P phases; its exponent is u = 0.6349. This surface intersects the 
surface JF hmiting the F phase at positive J^i as shown in Fig. 2. The model (|1|) 
exhibits the exact symmetry J7i —J\ [3- This symmetry is not respected in 
Fig. 2. However, we observe that simple block transformations would completely 
miss the F-AF transition in the 2D version of the model (P ||25|. On the surface 
between the F and the AF phases we find the fixed point (d); it has one relevant 
eigenvalue given at p = p* by A = with yj) = 2.72454. We interpret the point 
(D) as a discontinuity fixed point related to the F-AF first-order transition, which 
should be characterized by the value yr, = D = 3 |2^. If we maximize the free 
energy with respect to the discontinuity fixed point, we get yn = 1-78 at p = 0.31, 
which is the lowest value for which the discontinuity fixed point exists. The fact 
that this result is worse than the one obtained ai p = p* can be explained by saying 
that the LBRG transformation does not give good results when not ferromagnetic 
ground-states are involved. 

Numerical simulations of P, show the existence of a line of tricritical points 
on the F-P transition surface close to the J7i = plane; this line, at decreasing 
values of J74, ends in a Baxter point. Due to the limits of applicability of the LBRG 
transformation at small J7i, we cannot give reliable predictions on the structure 
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of the phase diagram in the region where the F-P and the AF-P surfaces meet. 
However, we have also studied the RG recursions on the hne L separating the 
domains of attraction of the fixed points (c) and (d) on the surface T . On the 
hne £, which is very close to the intersection of the AF-P with the JF surface, 
we find a fixed point (l) with two relevant eigenvalues, which annihilates with 
the discontinuity fixed point for j9 < 0.31. The free-energy of this fixed point is 
maximum at "p very close to p*, where the exponents of the two relevant eigenvalues 



are y\ = 1.59347 and t/2 = 0.09564 ||2^. The largest not relevant eigenvalue is 



A = 0.90395. The fixed point (l) is also reported in Table 1; it can be seen that it 
is very close to the fixed point (c). A realistic discussion of the phase diagram in 
the plane J^i = 0, where the F-P and the AF-P surfaces should meet, is given in 



28 . 



In Fig. 3 the phase diagram is shown in the particular case J'^ = J^2, which 
means jSc = in the surface representation. The paramagnetic phase, in accord 
with Monte Carlo results and differently from what comes out from mean field ap- 
proximation [^], extends at positive J'l towards zero temperature. This is related 
to the high degeneracy of the ground states in this region [^. 

A different representation of the phase diagram can be given in terms of the 
surface parameters Pl,Ps and Pc (see (|^)). In Fig.4 the F-P-AF transitions are 
shown in the plane Pl,Ps for different values of the curvature Pc- The F phase, 
at large values of Ps: describes configurations with diluted small surfaces. By 
decreasing the value of Ps, area is favoured to increase and, at the percolation 
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threshold, interfaces invade the system. However, it is still possible to distinguish 
between an inside volume wrapped up in interfaces and a different outside volume. 
By decreasing furtherly the value of Ps, if Pl is sufficiently large, at the Ising- 



like F-P transition, a random isotropic |29| phase is stable and the symmetry of 



the hamiltonian between inside and outside is restored. The AF phase can be 
intended as a droplet crystal. The limit (3l oo describes a gas of self-avoiding 
surfaces and is particularly relevant for physics. In Table 2 the critical values of 
Ps for self-avoiding surfaces are reported at different values of the curvature and 
compared with results from simulations. 

Open interfaces - model (||,|^. The gauge model @ at /5m = is dual to the 



3D Ising model |^. At small /3m it can be expanded on the dual lattice as an 
Ising model with many interactions. At the second order of this expansion the 
gauge model is mapped onto the model (|]). The LBRG transformation has been 
applied to study the F-P transition in the model (^. Then the results have been 
reported by formulas in the plane PmPg, as shown in Fig.5. At small Pm the 
critical line starts from the Pa axis at Pg = — | lntanh/5f^j^, with /5^^j^ = 0.23926. 
The continuous line in Fig.5 is the self-dual line, which is the line mapped onto 
itself by transformations @ , with respect to which the phase diagram has to be 
symmetric. Then the critical line at small Pm is mapped by eqs. (|) in the region 
at large Pg- The two lines meet at the point Pm = 0.241, /J^ = 0.719 on the 



self-dual line. In the phase diagram found by numerical simulations |T6|, the two 
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lines starting at = and at (3g = oo become first-order at tricritical points 
before meeting on tfie self-dual line. There, at a triple point, an other first-order 
line comes out towards greater values of /Sm on the self-dual line. This first-order 
line ends with a critical point at finite and positive values of (3g and /Sm- As 



discussed in an interesting aspect of the expansion (|]) is that the 4-spin 
interaction terms are expected to give tricritical points. However, a mean-field 
approximation of the model (^) gives tricritical points quite far beyond the 



triple point. Also our calculations suggest that the transition lines are continuous 
on the parts drawn in Fig.5, which correspond to the critical F-P transition in 
the Ising representation. Therefore, the relevance of the 4-spin interaction in (|^) 
is probably not sufficient to explain alone the existence of the tricritical points 
found in simulations. These results will be further commented in the next section. 



§4. Discussions and conclusions. 

We have applied the LBRG transformation to study the phase diagram of Ising 
models describing closed and open interfaces. Interacting closed interfaces can be 
naturally expressed as an Ising model, while open surfaces, originally written as 
a gauge model with statistical variables on the bonds, can be mapped on Ising 
models only in extreme regions of the phase diagram. At large the gauge 
model describes the interesting physical situation of almost-closed surfaces with 
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few defects. 

First consider the model (0) of closed interfaces. Results concerning the tran- 
sition on the nearest-neighbour axis are in good agreement with previous known 
results. Also the value of the Ising exponent v = 0.6289 is in excellent agreement 
with other numerical work. We expect that the critical surface has been found 
with a good approximation in the region close to the nearest-neighbour axis. Re- 
sults regarding the interesting case of self-avoiding surfaces have been reported in 
Table 2. 

Problems arise when the LBRG transformation is applied to study regions of 
the phase diagram where ordered not ferromagnetic configurations are relevant. 
In particular, the LBRG transformation does not take into account the J7i —J\ 
symmetry of the model (P which should give at low temperatures a first-order 
F-AF transition at jTi = 0. We find this first-order transition, but not at ^Ti = 
(see Figs. 2, 3). Moreover, our results cannot reliably describe the region close to 
the line where the F-P and the AF-P surfaces meet, which should be on the plane 
J7i = 0. However, for completeness, we have also given results concerning this 
region. 

The model (|I|) has been largely studied in D = 2 [Q, where RG transforma- 



tions taking correctly into account the ground-state structure have been consid- 



ered giving the expected topology of the phase diagram [£8|, |3ll]. We have tried to 
generalize the LBRG transformation in order to take into account the existence 
of antiferromagnetic ground states. Then we have considered a weight function 
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V{(j' ^ a) distinguishing between spins of different sublattices. For each cell marked 
by a cross in Fig.l the spin a' is coupled only to the four spins of one original 
sublattice (see eq. (|^) and Fig.l), in such a way that two nearest neighbouring 
spins a' are coupled to the spins a of different sublattices. Then all the interac- 
tion is moved into the dark grey cells of Fig.l and a RG transformation analogous 
to eq.(^) can be written in such a way to get a homogeneous hamiltonian with 
the same expression for any elementary cell. By this procedure we have obtained 
phase diagrams which exhibit the symmetry J7i —J\i but with a rather poor 
precision for the critical temperature on the nearest-neigbour axis and for the ex- 
ponent V. Moreover the tricritical points numerically found close to the plane 
iTi = are not obtained by this transformation. Therefore a complete RG study 
of the phase diagram of the model (|l|) in D = 3 is still an open question. 

In Fig. 5 we have presented the phase diagram of the self-dual 2^(2) gauge model 
found by applying the LBRG transformation to the model (|^). Our estimation 
of the critical lines is reliable especially in the region of validity of the expansion 
(|^), that is at small /?m and, by duality, at large close to the points where the 
model can be written as an Ising model with only nearest-neighbour interaction. 
Numerical simulations |16| predict that these lines become first-order before meet- 



ing on the self-dual line. By our methods, we cannot predict such a behaviour. 
Indeed, our results suggest that the transition line remain continuous for a long 
part beyond the self-dual line. Therefore, even if tricritical points could arise in 
model d^), that expansion is probably not useful to discuss the phase diagram of 
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the model (|]) close to the self-dual line, for which other methods are needed. In 
conclusions, provided all the discussed limitations, we can say that the application 
of the LBRG transformation to spin models describing lattice interfaces gives, in 
a relatively simple way, phase diagrams in many parameter spaces which are quite 
accurate especially in the region close to the nearest-neighbour axis. 
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Table Captions 

Table 1. 

The coordinates of the fixed point related to the F-AF-P transitions are reported 
for the value p = p* which maximizes the free energy of the critical F-P fixed 
point (C). The other symbols (f), (af), (p), (ac), (d) and (l) denote respectively 
the low-temperature ferromagnetic and antiferromagnetic fixed points, the high- 
temperature, the critical AF-P, the discontinuity fixed point between the AF and 
F phases, and the fixed point on the manifold separating the domains of attraction 
of the fixed points (c) and (d) (on the hypersurface limiting the F phase). The 
two squares on the left represent two parallel faces of an elementary cube of the 



lattice. The dots represent the spins taking part in a given interaction. In [19| the 
value of p maximizing the critical F-P fixed point free energy has been found to 
he p = 0.40343. This fixed point is symmetric (in the sense explained in the main 
text) and the 2-spins, 4-spins, 6-spins and 8-spins coordinates are respectively 
0.02097, 1.96 X 10-^ -7.69 x 10"^ and 2.15 x 10"^ 
Table 2. 

The critical values of j3s in the self-avoidness limit for different values of Pc- The 
Monte Carlo results are taken from Refs p and p. 
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Figure Captions 

Fig.l: A 2D representation of the LBRG transformation. The crosses indicate 
the spins cr'; the squares are the original spins a. The a-dependent terms in the 
hamihonian and in the weight function are moved into the grey squares. Full and 
empty squares represent spins of the two original sublattices. In the variant of the 
LBRG transformation described in §4 the a-dependent terms are moved into the 
dark grey cells of the lattice. 

Fig. 2: The phase diagram obtained by applying the LBRG transformation to 
the model (0). The horizontal surface at positive J2 separates the F phase at 
positive Ji from the AF phase. At lower values of J2 the F and the AF phases 
are separated by the P phase. 

Fig. 3: The phase diagram of the model (|l]) with J2 = Jj^. Fixed points of Table 
1 are also reported. 

Fig. 4: The phase diagram of the model ([^) in terms of the surface parameters 
(3s and jSi- The different curves refer from the right to the left respectively 
to the values of f3c = -0.2,-0.1,-0.04,0. The symbols F, AF and RI denote 
respectively the ferromagnetic, the antiferromagnetic and the random isotropic or 
paramagnetic phase. 

Fig. 5: The phase diagram of the self-dual Z{2) gauge model. The continuous 
line is the self-dual line. The dashed lines are critical lines found by applying the 
LBRG transformation to the model (|^). 
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Table 1 
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Table 2 



0.04 -0.04 -0.1 -0.2 

(3, LBRG 0.446 0.472 0.504 0.560 0.674 

Monte Carlo - O.SSSl^l O.SeoI^l 0.4701^1 0.5701^1 



